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ABSTRACT 

We examine nucleosynthesis in the electron capture supernovae of progenitor AGB stars with an O- 
Ne-Mg core (with the initial stellar mass of 8.8 Mq). Thermodynamic trajectories for the first 810 ms 
after core bounce are taken from a recent state-of-the-art hydrodynamic simulation. The presented 
nucleosynthesis results are characterized by a number of distinct features that are not shared with 
those of other supernovae from the collapse of stars with iron core (with initial stellar masses of more 
than 10 Mq). First is the small amount of 56 Ni (= 0.002 — 0.004 Mq) in the ejecta, which can be an 
explanation for observed properties of faint supernovae such as SNe 2008S and 1997D. In addition, the 
large Ni/Fe ratio is in reasonable agreement with the spectroscopic result of the Crab nebula (the relic 
of SN 1054). Second is the large production of 64 Zn, 70 Ge, light p-nuclei ( 74 Se, 78 Kr, 84 Sr, and 92 Mo), 
and in particular, 90 Zr, which originates from the low Y e (= 0.46 — 0.49, the number of electrons per 
nucleon) ejecta. We find, however, that only a 1 — 2% increase of the minimum Y e moderates the 
overproduction of 90 Zr. In contrast, the production of 64 Zn is fairly robust against a small variation 
of Y e . This provides the upper limit of the occurrence of this type of events to be about 30% of all 
core-collapse supernovae. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — supernovae: 
general — supernovae: individual (SN 1054, SN 1997D, SN 2008S) — nebulae: Crab 
Nebula 



1. INTRODUCTION 

Massive stars end their lives with core-collapse super- 
novae (SNe II/Ibc), which are the predominant sources 
of having enriched galaxies with the elements heavier 
than helium. The other type, thermonuclear supernovae 
(SNe la), also contributes to the enrichment of iron-peak 
elements, which is, however, absent in the early universe. 
The metals produced by core-collapse supernovae serve 
as diagnostic tools to uncover the chemical-enrichment 
history of the Galaxy from its poorly understood early 
stage to the present day. A reliable prediction of su- 
pernova yields has been, however, hampered by the yet 
unknown mechanism that causes the explosion. Pre- 
vious studies of supernova nucleosynthesis have relied 
upon a number of model parameters such as the explo- 
sion energy, the position that divides the ejecta and the 
remnant (mass cut), and the electron fraction (Y e , the 
number of electrons per nucleon). The production of 
each element in a supernova, in particular of those syn- 
thesized in the innermost ejecta, is sev erely affected by 
the choice of these parameter s (see, e.g.. lTominaga etal] 
120071 iHeger fc Woosfevll2008ft . It is obvious that nucle- 
osynthesis studies with self-consistently exploding mod- 
els are eventually needed to obtain reliable supernova 
yields. 

Recent one-dimensional simulations including ac- 
curate neutrino transport seem to exclude the pos- 
sibility of neutrino-driven explosions without the 

1 Institute for the Physics and Mathematics of the Universe, Uni- 
versity of Tokyo, Kashiwa, Chiba 277-8582, Japan 

2 Department of Astronomy, School of Science, University of 
Tokyo, Bunkyo-ku, Tokyo, 113-8654, Japan; wanajo@astron.s.u- 
tokyo.ac.jp, nomoto@astron.s.u-tokyo. ac.jp 

3 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild- 
Str. 1, D-85748 Garching, Germany; thj@mpa-garching.mpg.de, 
kitaura@mpa-garching.mpg.de 



2001 



Buras et al 



help o f multi-dimension al effec ts (iRampp fc Janka 
20001 iMezzacappa et all l200lt iLiebendorfer et al. 
[Thompson et alJ 12003b ISumivoshi et all 120051 : 
" 12006ft . at least for standard nuclear and 
neutrino physics input. An exception are the explosions 
of 8 — 10 Mq stars. A star in this mass range forms 
an electron-degenerate core consisting of oxygen, neon, 
and magnesium (O-Ne-Mg) during the final stage of 
its evolution (instead of an iron core in the case of 
more massive stars), and becomes a super asymptotic 
giant branch (AGB) star. Such a star ends its life 
either as an O-Ne-Mg white dwarf or a core-collapse 
supernova leaving behind a neutron star. For the latter, 
the collapse is induced by electron capture (what is 
called an "electron capture supernova" ) when the core 
mass grows to 1.38 Mq and the central density reaches 
4 x 10 9 g cm~ 3 . However, the uncertainties in mixing 
and mass loss during the evolution make it di fficult to 
draw a clea r line between these two channels (|Nomotol 
1984, 1987). Recent studies report that only the range 
close to the upper e nd (a mass r ange of < 1 Mq or ~ 4% 
of all supernovae, ISiessI [20071 iPoelarends etail [2008) 
leads to the explosion channel, although the range could 
be wider for lower metallicity stars. 

The structure of the O-Ne-Mg core is distinctively dif- 
ferent from the iron cores of more massive stars by the 
fact that it has a steep density gradient in the outermost 
layers, surrounded by an extrem ely extended, l oosely 
bound H/He envelope. Recently, iKitaura et al.1 (|2006l ) 
have obtained self-consistent explosions from the collapse 
of the O-Ne-Mg core in a stell ar progenitor w ith an initial 
mass of 8.8 M developed bv lNomotol (|1984D . Their one- 
dimensional simulations with a state-of-the-art, energy- 
dependent treatment of the neutrino transport are in 
fact the only recent models in the literature with suc- 
cessful supernova explosions in spherical symmetry for 
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standard nuclear an d weak interaction physics (see also 
iBurrows et alJl2~007l . for a similar result). The explosions 
are initiated by the rapid outward acceleration of the su- 
pernova shock when it encounters the steep density gradi- 
ent and fast decline of the mass accretion rate at the edge 
of the O-Ne-Mg core. They are powered by the neutrino- 
heating mechanism, which yields a low explosion en- 
ergy of ~l-2 x 10 50 erg (for a detailed discussion, see 
Uanka et alj |2008a). The new calculations are a revision 
of previous hydrodynamic results for t he same O-Ne-Mg 
core, namely, of prompt e xplosions ([Hillebrandt et all 
Il984t IWanajo et al.l 12003). powerful neutrino-driven 
explo sions (0.6 - 1.2 x 10 51 erg, iMavle fc Wilsonl 
1988) , and no explosions (IBurrows fc Lattimerl 119851 : 
Baron. Cooperstein. fc Kahanalll987j T 

The purpose of this paper is to present nucleosynthesis 
results for the first 810 ms of the neutrino-driven explo- 
sion of a collapsing O-Ne-Mg core (electron capture su- 
pe rnova), using the ther modynamic trajectories obtained 
bv iKitaura et al.l |2006). Employing a self-consistently 
calculated explosion model with sophisticated neutrino 
transport is of particular importance not only for the nu- 
cleosynthesis study itself, but also for a couple of other 
aspects. On the one hand, the mass range of 8-10 M© ac- 
counts for about 30% of all the core-collapse supernova 
events, if all the range leads to the explosion channel. 
The electron capture supernovae can thus be potentially 
significant contributors to the Galactic chemical evolu- 
tion of some species. On the other hand, explosions from 
these progenitors have been proposed as a possible expla- 
nation for the infe rred low explosion energy of SN 1054 
(Crab supernov a. iNomoto et al.l |1982t I Davidson et al.1 
I1982Q (see also iHillebrandtl I1982L for the 10 M pro- 
genitor with an iron core) as well as for the small 56 Ni 
amount estimate d for some low-luminosity supernovae 
(e.g.. SN 1997D.[Turatto et a"Dll998t IChugai /fc Utrobml 
I2000l: iBenetti et al. 2001; H endry et al.l 120051 ) . A newly 
identified class of luminous transients like SN 2008S, 
whose progenitors were deeply dust-enshrouded mas- 
sive stars, is also suggested to be electron captur e su- 
per novae dPrieto et al.1120081 : iThompson et al.1120081 ) (see 
also lPastorello et al.ll2007f) . 

Recently. iHoffman et al.1 (|2008l ) have investigated the 
nucleosynthesis in elect ron capture superno vae, using an 
explosion calculation of Uanka et al.l (l2008al) that i s very 
similar to the explosion models of IKitaura et all (|2006f ) 
(but was computed wit h slightly different input physics, 
see Uanka et al.ll2008 a). Their study was aimed at de- 
termining whether the conditions are suitable for an r- 
process. They found no r-process formation and a severe 
overproduction of 90 Zr, which is also seen in our results. 
In our work, however, we perform nucleosynthesis calcu- 
lations in much more detail, taking into account a num- 
ber of possible uncertainties. In the follow i ng se ction 
(§ 2) the explosion models of IKitaura et al.l (|2006l ) and 
the methods for the nucleosynthesis calculations will be 
described. The nucleosynthesis results for the original 
model and for various modifications of it will be pre- 
sented in § 3. We will discuss the question whether elec- 
tron capture supernovae can be significant contributors 
to Galactic chemical evolution in § 4 and will address the 
possibility that they are the origin of some low-energy su- 
pernovae in § 5. Finally, our summary and conclusions 
will follow in § 6. 



2. EXPLOSION MODEL AND REACTION NETWORK 

IKitaura et alj (|2006f ) have simulated the co llapse of 
the O -Ne-Mg c ore (with the mass of 1.38 M , iNomotol 
QUI) (see also iMivaii et al.lfl980| ). employing two dif- 
ferent ("soft" and "stiff") nuclear equations of state 
(EoSs ). In this study, we adopt the result with the softer 
EoS (jLattimer fe Swestvl ll99lL LS) as the "standard" 
model (labelled ST), which is qualitatively similar to that 
with the stiffer one (W H; Wolff & Hillebrandt EoS in 
iHillebrandt et al.lll984D . The model with the WH EoS is 
used later for the purpose of comparison (§ 3.2). Figure 1 
shows the evolution of the radius, density, temperature, 
and Y e for selected mass elements in model ST as a func- 
tion of post-bounce time t p b- The ejecta in model ST are 
split into 29 mass shells with the first ejected zone having 
a mass of about 4 x 10~ 4 M© and each of the other ones 
of about 5 x 10~ 4 M© (1.39 x 10" 2 M Q in total). 

The nucleosynthesis yields in each mass shell are ob- 
tained in a post-processing step by solving an extensive 
nuclear reaction network code. The network consists 
of 6300 species between the proton- and neutron-drip 
lines predicte d by the recent fully microscopic mass for- 
mula (HFB-9. iGorielv et al.ll2005l ), all the way from sin- 
gle neutrons and protons up to the Z = 110 isotopes. 
All relevant reactions, i.e. (n,j), (p, 7), (a, 7), {p,n), 
(a,n), (a,p), and their inverses are included. The ex- 
perimental data, whenever available, and the theoreti- 
cal predictions for light nuclei (Z < 10) are taken from 
the REACLIB 4 compilation. All the other reaction rates 
are taken from the H auser-Feshbach rates of BRUSLIB 5 
( Aikawa et al.ll2005[) making use o f experimental masses 
([Audi. Wapstra. & Thib ault 200 3]) whenever avai lable or 
the HFB-9 mass predictions (jGorielv et al1l2005f ) other- 
wise. The photodisintegration rates are deduced from 
the reverse rates applying the reciprocity theorem with 
the nuclear masses considered. The weak and interme- 
diate screening c orrections to cha r ged p article reactions 
are adopted from iGraboske et al.1 (|1973l ) . 

The /3-decay rates a re taken from the gross theory 
predi ctions (GT2, iTachibana. Yamada. fc Yoshidal 
I1990D obtained with the HFB-9 predictions (T. 
Tachibana 2005, private communication). Elec- 
tron capture reac tions on free nucleons an d on 
heavy nuclei dFuller. Fowler, fc Newman! 119821 : 
lLanganke fc Martinez-Pinedd l2001h as well as rates 
for neutrino capture on free nucleons and 4 He and 
for neutrino s p allatio n of free nucleons from 4 He 
(IWooslev et al.l 119901: iMcLaughlin. Fuller, fc Wilsonl 
I1996D are also included. In contrast, neutrino-induced 
reactions of heavy nuclei are not taken into account in 
this stu dy, but they are expected to make only minor 
effects (jMever. McLaughlin, fc Fullerl 119981 ). Figure 2 
shows the luminosities and mean energies for neutrinos 
of all types as functions of the post-boun ce time t n h in 
model ST; these results are taken from IKitaura et al.l 
(2006) and used for calculating the rates of neutrino- 
induced reactions. More preciously, we should apply 
those quantities in the co-moving frame of the fluid 
with corrections for the gravitational redshift and for 
Dopplcr shifting due to fluid motion. We neglect such 
effects, which are important on the one hand only 

4 http:/ /nucastro.org/reaclib.html. 

5 http:/ /www. astro. ulb.ac.be/Html/bruslib. html. 
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Fig. 1. — Radius (a), density (b), temperature (c), and Y e (d) 
as a function of post-bounce time for material ejected from the 
collapsing O-Ne-Mg core in model ST. The trajectories are colored 
in green, blue, and red for Y e < 0.470, 0.470 < Y e < 0.500, and 
Y e > 0.500 (values at the end of simulation), respectively. 

when the fluid is relatively close to the neutrino sphere 
(< several 10 km), where the temperature is still higher 
than Tg = 9 (see below). On the other hand, at large 
distances, where the expansion velocities of the gas 
are larger, neutrino interactions become essentially 
irrelevant. 

Each nucleosynthesis calculation is initiated when the 
temperature decreases to Tg = 9 (where Tg = T/10 9 K). 
In the first ejected trajectory, the highest temperature is 
Tg w 7 (Fig. 1), which is taken to be the initial condition 
for this case only. At such high temperatures, the com- 
position is in the nuclear statistical equilibrium (mostly 
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Fig. 2. — Neutrino luminosities (top) and mean neutrino energies 
(bottom) as functions of post-bounce time for electron (dotted line), 
anti-electron (solid line), and heavy-lepton (dashed line) neutrinos. 
The data are given for an observer at rest at 400 km from the center. 
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Fig. 3. — Initial electron fraction Y £il for material ejected from 
the core as a function of the enclosed mass, M r . Different colors 
correspond to models ST, FP3, FM3, MX, and WH as denoted in 
the panel (see text). 

free nucleons and few a particles), which is realized im- 
mediately after the calculation starts. The initial com- 
positions is then given by X n = 1 — Y e \ and X p = Y e i , 
respectively, where X n and X p are the mass fractions of 
free neutrons and protons, and Y e ^ is the initial electron 
fraction at Tg = 9 (Fig. 3, black line for model ST). 

3. NUCLEOSYNTHESIS RESULTS 

In subsection 3.1, we will present the n ucleosynthesis 
result s for the unmodified model (ST) of iKitaura et al.l 
(2006), and the corresponding information for variations 
of model ST will then be discussed in the following sub- 
sections. In each model, the nucleosynthetic yields for all 
the trajectories are mass-integrated over the ejecta-mass 
range. 

3.1. Unmodified Model 
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Fig. 4. — Mass fractions of isotopes ( after decay) in the ejecta of 
model ST relative to their solar values (Loddcrs 2003) as functions 
of mass number. The abundances smaller than X/Xq < 40 are 
omitted. The even-Z and odd-Z isotopes are denoted by open 
circles and triangles, respectively. The p-nuclei are represented 
with filled symbols. The solid lines connect isotopes of a given 
element. The dotted horizontal lines indicate a "normalization 
band" between the largest production factor and a factor of ten 
smaller than that, along with the median value (dashed line). 
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Fig. 5. — Mass fractions of elements (after decay) in the ejecta of 
model ST relative to their solar values (Loddcrs 2003) as a function 
of atomic number. The even-Z and odd-Z elements are denoted 
by open circles and triangles, respectively. 

The nucleosynthesis results of the unmodified model 
(ST) are shown in Figures 4 (isotopes) and 5 (elements). 
Both plots present the overproduction factors defined by 
the mass f ractions in the ejecta with respect to their so- 
lar values (|Lodderdl2003D . The even-Z and odd-Z species 
are denoted by circles and triangles, respectively. In Fig- 
ure 4, the isotopes for a given element are connected by 
lines, and the abundances smaller than X/Xq < 40 are 
omitted. The dotted horizontal lines indicate a "nor- 
malization band" between the largest production factor 
( 90 Zr and Zr in Figs. 4 and 5, respectively) and a factor 
of ten smaller than that, along with the median value 
(dashed line). This band is taken to be representative 
of the uncertainty in the nuclear data involved. In the 
following, we consider that electron capture supernovae 
can be contributors to the solar (or Galactic) inventories 
of the species located within the normalization band. 

Figure 4 indicates that, along with the marginal ones, 
model ST can account for the production of 64,66 Zn, 
70 Ge, 74 < 76 Se, 78 > 80 Kr, 84 Sr, 90 Zr, and 92 Mo, where all the 
light p-nuclei (filled circles) up to A = 92 are included. 
However, only a few elements (Zn, Ge, Y, and Zr) fall 
into the normalization band (Fig. 5), since the p-nuclei 
comprise only small fractions of a given element (0.89%, 
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Fig. 6. — Mass fractions of several important isotopes (after 
decay) in the ejecta of model ST (top), and these mass fractions 
relative to their solar values (bottom) as functions of the enclosed 
mass M r . The vertical dotted lines indicate the mass coordinates 
where the initial Y B>i is 0.470, 0.480, and 0.490 (see Fig. 3). 

0.35%, 0.56%, and 14.8% for 74 Se, 78 Kr, 84 Sr, and 92 Mo, 
respectively). Our resul t for model ST is in r easonable 
agreement with that in iHoffman et alj {2008, Fig. I) 6 , 
and some differences are due to the slightly higher mini- 
mum Y e in our model (§ 3.5). 

Figure 6 shows the mass fractions of some important 
isotopes (after decay) for model ST (top) and these mass 
fractions relative to their solar values (bottom) as func- 
tions of the enclosed mass M r . In each mass shell, 
the dominant heavy isotope is either 56 Fe (produced in 
the form of 56 Ni) or 58 < 60 Ni, depending on F eji . In the 
neutron-rich ejecta, the nuclear products in the single 
quasi-statistical-equilibrium (QSE) cluster have peaks at 
A ~ 60 and A _~_90 because of the strong binding at 
N = 28 and 50 (|Mever. Krishnan. fe Clavtoii i998). The 
abundance of 90 Zr is maximal in the mass shell with 
the lowest Y e \ = 0.464. In addition, the matter with 
Y e< i = 0.46 — 0.49 is so proton-rich (compared to the 
average of the stable isotopes in the vicinity, e.g., the 
proton fraction of 90 Zr is 0.44) that the successive pro- 
ton captures lead to the production of light p-nuclei up 
to A = 92. T his can also be seen i n the early ne utrino- 
driven winds (jHoffman et al.lll99d IWanajdl2006t ). 

It should be noted that such a nucleosynthesis result 
that stems from the low-Y e (~ 0.47 — 0.49) matter in the 
early ejecta may be a unique characteristics of collapsing 
O-Ne-Mg cores. In the present model, the explosion sets 
in immediately after core bounce and the ejecta rapidly 
expand, where the ejection of low Y e matter seems un- 

6 The total ejecta ma ss including th e oute r H/He envelope is 
taken to be 1.263 M in IHoffman et all |[200ST ). while we consider 
only the calculated zones with 1.39 X 10 -2 M Q . This leads to the 
hundred times larger values of X/Xq in this paper. 
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avoidable (see Fig. 1). For more massive progenitors 
(> 10 Mq), the explosion is expected to be more delayed 
and to eject neutrino-processed m atter, where the bulk 
of the ejecta may have Y e > 0.5 (|Frohlich et al.l [20061: 
iBuras et alJl2006h . 

As will be discussed in § 4, the maximum overproduc- 
tion factor of X/Xq = 5.7 x 10 5 for model ST (Fig. 4) 
poses a severe constraint on the occurrence of this type 
of events to be no more than 1% of all core-collapse su- 
pernovae. This is in agreement with the conclusions by 
Hoff man et al.l (120081). In the sim ulation of lKitaura et al.l 
(|2006h (also iJanka et afll2008ah . the deceleration of the 
shock in the outer envelope slows the expansion of the 
ejecta only slightly. Therefore we do not expect any sub- 
stantial fallback of the once ejected matter onto the rem- 
nant, as it is presumed to take place in the case of more 
massive progenitors (e.g., lUmeda fc Nomotol I2002D . In 
the following subsections we thus explore possible mod- 
ifications to model ST, which instead of fallback might 
provide a solution for moderating the extremely large 
overproduction of, in particular, 90 Zr. 

3.2. Uncertainty in the Equation of States 

First, we examine nucleosynthesis in the explosion 
model calculated with the WH EoS in iKitaura et al.l 
(2006). This model (hereafter WH) has thermody- 
namic trajectories very similar to those of model ST, 
but a slightly larger ejecta mass (1.64 x 10~ 2 M Q ) and 
a s lightly lower mini mal Y e ^ (= 0.462). Another EoS 
by IS hen et al.1 (|1998l ), which is also currently available 
for core-collapse simulations, falls between the ST and 
WH EoSs in terms of its stiffness and a variety of re- 
sults of core-collapse simulations, e.g. the radii of shock 
formation a nd stagnation or the size of the neutrino lumi- 
nosities (see lJanka et al.ll2008b| ). Therefore, we suspect 
that a comparison of the nucleosynthesis results between 
models ST and WH well brackets the uncertainties aris- 
ing from different EoSs. 

All nucleosynthesis calculations are repeated with the 
thermodynamic trajectories and the initial compositions 
deduced from the Y^i-mass profile (Fig. 3) of model 
WH. Figure 7a shows the mass fractions (bottom) in the 
ejecta for models ST and WH, along with their ratios 
Awh/Ast (top), as functions of the mass number A. 
Differences exceeding a factor of 2 (indicated by dotted 
lines) can be seen for the light species with A < 20. 
However, the differences are well below a factor of 2 for 
the dominant species in the vicinity of the QSE peaks 
at A 4,60, and 90. We therefore conclude that the 
uncertainties arising from different nuclear EoSs are not 
of great importance, at least for the currently available 
versions of EoSs. We note, however, that the WH EoS 
leads to a 60% larger 56 Ni mass than that for the LS 
EoS (Table 1). This suggests that EoSs play a crucial 
role to precisely determine the 56 Ni ejecta mass from an 
electron capture supernova. 

3.3. Uncertainty in the Nuclear Reaction Rates 

In all the present calculations, the nucleosynthetic 
flows proceed along or in the vicinity of the /3-stability 
line, forming the QSE peaks at A ~ 4, 60, and 90. The 
nuclear masses of relevance, which are the most impor- 
tant facto rs in this case, are measured wit h high enough 
accuracy (|Audi. Wapstra. fc Thibault]|2003D . Therefore, 
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Fig. 7. — Mass fractions (bottom of each panel) in the ejecta 
for models ST (blue) and (a) WH, (b) RT, (c) nou (red; see the 
text), and their ratios (top of each panel), as functions of the mass 
number. In the top panels, a factor of 2 difference is indicated by 
dotted lines. 

we do not expect a sizable change of the nucleosynthesis 
result arising from the uncertainties in nuclear reaction 
rates. 

As a test, we repeat all the nucleosynthes is calculations 
of m odel ST with the theoret ical rates (lAikawa et al. 
l2005f) replaced by those of iRauscher fc Thielemann 
(2000). Both data sets are based on the statisti- 
cal (Hauser-Feshbach) approach with the experimental 
masses whenever available, but with different theoretical 
masses and different nuclear level densities. The results 
(labelled by RT) are compared to those of model ST in 
Figure 7b. Factors of 2 — 3 difference can be seen for 
the nuclei between A = 27 and 45 (Fig. 8, top), which 
have, however, very small mass fractions (bottom). The 
differences for dominant species around A « 4, 60, and 
90 are well below a factor of 2. 

3.4. Effects of Convection 
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A potentially important effec t that is lacking in the 
one-dimensional simulations of [Kitaura et alJ (|2006h is 
convective mixing of the ejecta (see the results of a two- 
dimensional simulation °f a collapsing O-Ne-Mg core in 
Uanka et alJl2008bf ). Fi gure 8 shows the profiles of tem- 
perature (top), entropy (middle), and Y e (bottom) versus 
mass at early times (i p b = 126, 150, 176, 200, and 228 ms) 
and at the end of the simulation (t p b = 810 ms). We find 
a negative entropy gradient forming at a temperature 
high enough for a-processing (Tg > 3, which is indicated 
by the dashed line in the top panel of Fig. 9). This can 
cause convective overturn and one might speculate that 
this could moderate the neutron-richness to some extent 
before the freezeout. 

As a limiting case, we present the result with the initial 
composition being determined by the Yg^-mass profile 
shown by the green line in Figure 3 (hereafter, model 
MX). We assume here that the ejecta between 1.367 and 
1.375 Mq (indicated by vertical dotted lines in Fig. 9) 
get completely mixed on microscopic scales to have a 
mass-averaged constant Y e .i (— 0.480). We find that the 
species having the largest overproduction, i.e., 90 Zr in 
model ST, is replaced by 64 Zn with ten times smaller 
value in model MX (Fig. 9). The largest overproduction 
of Zr (as element) is replaced by that of Zn (Fig. 10). 
Because of the increased minimum value of Y" e ,i (Ye, mm = 
0.480) the overproduction of 90 Zr becomes unimportant. 

It should be noted, however, that a corresponding two- 
dimensional simul ation (carried out u ntil 262 ms after 
core bounce; see iJanka et al.1 (|2008bl ) and Miiller and 
Janka 2008, in preparation) with the same initial O- 
Ne-Mg core does not show any such small-scale mixing. 
Instead, the accreted post-shock gas makes one quick 
overturn (fully developed convection with many over- 
turns does not occur), after which the rising material in 
Rayleigh- Taylor mushrooms is directly ejected and self- 
similar expansion is quickly established. Hence, there 
may not be sufficient time for mixing and homogeniza- 
tion on small scales. As will be discussed in § 3.5, how- 
ever, a slight increase of Ye.min ( e -g-j AY e ,i = 0.004; from 
0.464 to 0.468) is already enough to significantly moder- 
ate the overproduction of 90 Zr. Therefore, partial mixing 
of the ejecta induced by convection might be sufficient to 
cure this overproduction problem. 

The negative entropy gradient may have a variety 
of other effects. One is that the convective overturn 
could stretch the mean duration of the neutrino irra- 
diation of ejected matter and could also lead to re- 
gions of ejecta with higher Y e \. In fact, a recent two- 
dimensi onal simulation wit h a more massive progenitor 
(15 Mq, iBuras et all 120061 ) shows that the bulk of the 
ejecta turns out to be proton-rich (Y e > 0.5). In the 
present case, however, the convective overturn may not 
have exactly the same effects as in more massive pro- 
genitors, because the steep density gradient of the O-Nc- 
Mg core with its transition to an extremely dilute H-rich 
envelope makes the core structure distinctively different 
from that of more massive stars. As a consequence, the 
supernova ejecta accelerate much faster than in more 
massive stars, and the convective pattern in the over- 
turn region freezes out in the self-similar expansion more 
quickly. Therefore convective mixing is rather inefficient 
and, moreover, it is possible that the neutrino-heated 
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Fig. 8. — Profiles of temperature (upper), entropy (middle), 
and Y e (bottom) versus enclosed mass M r at early times (t pb = 
100 — 230 ms) and at the end of simulation (tpb = 810 ms). See the 
text for the meaning of the vertical dotted lines and the horizontal 
dashed line. 

bubbles rise so rapidly away from the neutrino-sphere 
that clumps of low- Ye material get ejected, a possibility 
that may be anticipated from the bottom panel of Fig- 
ure 8 (e.g., the red and green lines there). This might in- 
crease the overproduction of other N = 50 species ( 86 Kr, 
87 Rb, 88 Sr, 89 Y). In fact, some neutron-rich blobs (down 
to Y e ~ 0.41, but with a tiny mass) are found in the pre- 
liminary results of the mentioned two-dimensional simu- 
lation (Miiller and Janka 2008, in preparation, see also 
IJanka et al.ll2008bh . 

It should be noted that even two-dimensional models 
with their constraint of axisymmetry may not yield a 
sufficiently accurate mass distribution of the ejecta as a 
function of Y e . Such a sensitive information might ul- 
timately require three-dimensional simulations to allow 
for reliable conclusions concerning nucleosynthesis yields 
and production factors. 

3.5. Small Y e Variation 

From Figure 6 we conclude that the large overproduc- 
tion of 90 Zr in model ST is mainly due to the neutron- 



rich ejecta with Y ej i between Y ej 



0.464 and 0.470. 



The demonstration in § 3.4 (Fig. 9) shows that boost- 



Nucleosynthesis in Electron Capture Supernovae 



7 




mass number 



Fig. 9. — Same as Figure 4, but for model MX. 
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Fig. 10. — Same as Figure 5, but for model MX. 

ing Y^min to 0.480 indeed removes the overproduction of 
90 Zr. Motivated by this result, we explore in this sub- 
section that how much increase of Ye. m i n could cure this 
extreme overproduction problem. We repeat that the 
y ej i values are obtained from state-of-the-art simulations 
wi th sophisticated, ener gy-dependent neutrino transport 
in iKitaura et al.l (|2006[ ). Nevertheless it is a great chal- 
lenge to determine Y e to an accuracy of a few percent, 
because the neutron-to-proton ratio in the ejecta is estab- 
lished by a delicate competition of electron neutrino and 
antineutrino captures and their inverse reactions, in par- 
ticular around the radius where the kinetic equilibrium 
of these processes breaks down because the rates become 
slower than the expansion rat e of the accelerating ejecta 
(for a detailed discussion, see lFrohlich et al.1 12006). Be- 
sides the limited numerical resolution (in particular of 
the neutrino energy spectra) a variety of other effects 
can be imagined to imply uncertainties at the percent 
level, for example the potential effects of convective mix- 
ing on microscopic scales (§ 3.4), future refinements in 
the employed microphysics (e.g., EoSs, neutrino interac- 
tion rates, electron captures rates) , possible effects due to 
nonstandard neutrino properties (e.g., flavor oscillations 
in the supernova core), and the uncertainties associated 
with the stellar evolution calculations of 8 — 10 Mq stars. 

In order to test small variations of i^mim we examine 
nucleosynthesis for the trajectories of model ST but with 
the original Y e ^ profile replaced by Y^i+(0.50Q— Y e ;\) x /• 
The multiplicative factor / is taken to be 0.1, 0.2, and 
0.3 (hereafter, models FP1, FP2, and FP3). The Y e;i - 
M r profiles of these models are similar, but each model 
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Fig. 11. — Same as Figure 4, but for models FP1 (top), FP2 
(middle), and FP3 (bottom). 

has a different Ye, mm ( see Fig. 3 for model FP3). The 
values of Y eymin for these models are slightly increased 
compared to that of model ST to be 0.468, 0.471, and 
0.475, respectively (Table 1). 

We find that the production of 90 Zr is extremely sen- 
sitive to Ye.min (Figure 11). The increase of Y emin by 
only 1 - 2%'(AF e ,i = 0.004 and 0.011 in models FP1 and 
FP3, respectively) reduces the production factor of 90 Zr 
by roughly two orders of magnitude. Other species like 
64 Zn, 74 Se, and 78 Kr then possess the largest production 
factors. In the case of model FP3 these have one tenth of 
the values in model ST. Figure 12 indicates that electron 
capture supernovae could be the dominant sources of the 
elemental abundance of Zn, if such a slightly higher Y^; 
was correct. 

The results for / = —0.1, —0.2, and —0.3 (hereafter, 
models FM1, FM2, and FM3) are shown in Figures 13 
and 14. The values of Ye,min are then decreased to 0.460, 
0.457, and 0.453, respectively (Table 1). The Y e<i -M r 
profile for model FM3 is also displayed in Figure 3. The 
overproduction of 90 Zr becomes more serious in these 
models, and other N — 50 species ( 88 Sr and 89 Y) en- 
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Fig. 12. — Same as Figure 5, but for models FP1 (top), FP2 
(middle), and FP3 (bottom). 

ter into the normaliza tion band in t he ca se of model 
FM3. The results of IHoffman et all (|2008t Fig. 1) re- 
semble those for our model FM3 rather than those for 
model ST. Th i s is p robably explained by the fact that 
IHoffman et al.l (|2008f ) perform their evaluation for ejecta 
with a value of Y^min = 0.454, which is close to Y eim ; n in 
model FM3 (= (1453), but slightly smaller than V e , m in 
in model ST. The small dif f erence between the ejecta 
conditions in IHoffman et al.l (|2008l ) and those of model 
ST originates from a different density structure assumed 
for the dilute H/He envelop e around the c o llapsi ng O- 
Ne-Mg core consi dered by IHoffman et al.l (|2008l ) (see 
iJanka et al.l 120 08a. for more detail). 

The results presented in this subsection imply that the 
ov erproduction of 90 Zr m ay not be as serious as reported 
bv lHoffman et alJ (|2008[ L because an increase of Y einnn by 
only 2% cures this problem. This is a significant improve- 
m ent compared to the situ ation in the older simulations 
bv lMavle fc Wilson! lfl988h . where Y e m in was found to be 
around 0.4 and a substantial change (~ 20%) of Ye,min 
was necessary to avoid the overproduction of N = 50 
nuclei. 
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Fig. 13. — Same as Figure 4, but for models FM1 (top), FM2 
(middle), and FM3 (bottom). 

3.6. vp-process 

The j/p-process (or neutrino- induced rp-process) is now 
believed to be a promising nuclear process that synthe- 
sizes light p-nuclei up to A ~ 1 10 (|Frohlich et al.l 120061 : 
iPruet et al.ll2006UWanaidl2006f ). In this process, a frac- 
tion of free protons are converted to neutrons by neutrino 
capture in the early proton-rich supernova ejecta. The (3- 
waiting points on the classical rp-process path (e.g., 64 Ge 
with the half- life of 1.06 minutes) are then bypassed via 
much faster neutron capture. The current models con- 
tain proton-rich ejecta (up to 0.53, Figs. 1 and 3), where 
one may expect the occurrence of the i^p-process. 

The ^p-process plays, however, no role in producing 
the p-nuclei, although all the examined models in this 
study include neutrino- induced reactions on free nucleons 
and a particles. In order to test the effect of neutrinos, 
the result without neutrino-induced reactions (labelled 
not/, otherwise using the same input as in model ST) 
is compared with that of model ST (Fig. 7c). We find 
that neutron capture, which is absent without neutrino- 
induced reactions, diminishes the nuclei with A < 20. 
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Fig. 14. — Same as Figure 5, but for models FM1 {top), FM2 
(middle), and FM3 (bottom). 

The mass fractions of these isotopes are, however, very 
small compared to the dominant species at A rj 4, 60, 
and 90. No substantial differences can be seen for the 
nuclei with A > 50 because of the moderate proton- 
richness (up to Yg.i = 0.53), moderate entropy (up to 
30 /cb per nucleon, where fee is the Boltzmann constant), 
and the fast expansion of the ejecta. The proton-rich 
ejecta quickly cool down below Tg « 2 — 3, which is the 
relevant temperature range for the ^p-process to take 
place (i.e., the proton captures are fast and their inverse 
reactions slow) . The i/p-process might be efficient only in 
core-collapse supernovae from more massive progenitors 
(e.g., > 15 Mq), which have a shallower density gradient 
at the core edge and a denser envelope and thus their 
explosions develop in a different way with different con- 
ditions for nucleosynthesis. 

4. CONTRIBUTION TO GALACTIC CHEMICAL 
EVOLUTION 

We now discuss a possible contribution of electron cap- 
ture supernovae to Galactic chemical evolution. First, 
we suppose that model ST, which has the largest over- 



production factor for 90 Zr (Fig. 4), is representative 
of this type of supernovae. Let us assume that elec- 
tron capture supernovae produce all 90 Zr in nature, and 
the other supernovae from the progenitors more mas- 
sive than 10 Mq produce 16 with a typical amount of 
Moth CI ( 16 0) = 1.5 Mp per event. Here, M othcr ( 16 0) 
is taken to be the initi al-mass-function averaged yield 
of iNomoto et all (|2006f) between 13 and 40 M Q (solar 
metallicity models) . The contribution of electron capture 
supernovae to the Galactic 16 is negligible (Table 2). If 
we assume the number fraction of electron capture su- 
pernovae relative to all core-collapse supernovae to be 
/*, we have the relation 

/* _ *( 9O Zr) /X(i6Q) 
1 - /. ~ M( 9 °Zr)/M other ( 16 0) ~ U - U ^' W 

where X( 16 Q) ^ = 6.6 x 10" 3 , AT( 90 Zr) Q = 1.5 x 10~ 8 
(lLoddersll2003h . and M( 90 Zr) = 1.2 x 10~ 4 M© (Table 1). 
Therefore we expect the frequency of electron capture 
supernovae to be no more than 1% of all core-collapse 
events, when taking into account that there are also other 
sourc es of 90 Zr (81% is from the s-process. iBurris et al.l 
2000). As discussed in § 3.5, however, the production of 
90 Zr is extremely sensitive to Y e< i, changing more than 2 
orders of magnitude with only 2% variation of Y Bj \ (Ta- 
ble 1). Thus, we do not consider the overproduction of 
90 Zr to give a tight constraint on the occurrence of elec- 
tron capture supernovae. 

Instead, we propose the abundance of 64 Zn, whose 
production is insensitive to small variations of Y e ^ (Ta- 
ble 1), to serve as a strong constraint on the occurrence 
of this type of supernovae. As an example, we con- 
sider model FP3 (Fig. 11; bottom) as representative of 
electron capture supernovae, in which the largest over- 
production is shared by 64 Zn, 74 Se, and 78 Kr. The ori- 
gins of these isotopes have not been well identified, al- 
though the latter two p-isoto pes might be prod uced to 
some ex tent by the 7-proces s (iRavet et aTl 11995) or t^p- 
proce ss (jFrohlich et alj|2006t iPruet et al.l 120061 : IWanaid 
2006) in core-collapse supernovae. We therefore can as- 
sume that all 64 Zn (and 74 Se, 78 Kr) in nature is pro- 
duced by electron capture supernovae. In this case, equa- 
tion (1) with 90 Zr replace d by 64 Zn gives /* = 0.28, where 
AT( 64 Zn) = 1.1 x 10~ 6 (ILoddersll2003h and M( 64 Zn) = 
6.5 x 1O~ 4 M (Table 1). We thus conclude that the up- 
per limit of the frequency of electron capture supernovae 
is around 30% of all the core-collapse events. This is in 
good agreement with the upper limit est imated from a re- 
cent study of stellar evolution (~ 20%. iPoelarends et al.1 
I2008T) as well as the estima te for the rate of SN 2 008S- 
like transients (~ 20 - 30%. [Thompson et~all [2008) . Our 
result implies that a significant fraction of 8 — 10 Mq 
is allowed to enter into the supernova channel from the 
nucleosynthetic point of view. 

For a practical use, the yields of all stable (Table 2) 
and some unstable (Table 3) isotopes for models ST and 
FP3 are presented. For the study of Galactic chemi- 
cal evolution, one can use those of model FP3 as rep- 
resentative of electron capture supernovae, keeping in 
mind that the Y e ^ profile (i.e., the initial composition) 
of this model was slightly modified. Nevertheless, it 
would be interesting to see if this type of supernovae 
was the dominant source of Zn and of the most myste- 
rious p-nucleus 92 Mo. In fact, all the previous models 
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of supernova nucleosynthesis ( with typical explosion en - 
ergy), except for hypernovae (|Umeda fc Nomotoi r2002). 
have failed to explain the solar inventory of 64 Zn (which 
is the most abundant isotope of this element). We do 
not expect a substantial contribution to the abundance 
of Zn from the later (t v h > 1 s) n eutrino-driven wind 
(Hoff man et"all Il996t IWana id 120061. whose total mass 
is no more t han 0.001 Mq with 4 He being the dominant 
species (e.g., lWanaidl2007l ). The p-isotope 92 Mo can be 
synthesized by the up- process to some extent, but still 
falls a few times short of what is expected fro m the neigh- 
boring jp-nuclei ( 84 Sr, 94 Mo, and 96,98 Ru, iPruet et al.l 
l20l^lWanai^l200^ . 

5. ORIGIN OF FAINT SUPERNOVAE? 

Our result confirms that electron capture supernovae 
produce very little 56 Ni (rj 0.002 - 0.004 Mq, Table 1) 
compared to ~ 0.1 Mq in the cas e of more massive pro- 
genitors (e.g.. lNomoto et al.ll2006l ). This is a consequence 
of the small ejecta mass (= 0.0139 M Q without the H-rich 
envelope) and also of the neutron-richness of the bulk of 
the ejecta (Y e j < 0.49 for mass shells in the range of 
1.367 - 1.373 M Q , Fig. 6), because of which 56 Ni is not 
a dominant species to be produced. The yield of 56 Ni, 
which is mainly created in matter with Y e \ ~ 0.5, is in- 
sensitive to a small variation of Y e (Table 1). We also 
do not expect a significant contribution from the later 
(i p b > 1 s) neutrino wind. The uncertainty arising from 
the nuclear EoS seems larger than that from Y e , being 
about a factor of 2 (Table 1) for the currently available 
versions of EoSs. The uncertainty might become larger if 
a variety of EoSs are available in the future. The convec- 
tive motion near the mass cut may also affect the 56 Ni 
mass (§ 3.4), which needs multi-dimensional simulations. 

The expected small amount of 56 Ni as well as the 
low explosion energy of electron capture supernovae have 
been proposed as an explanation of t he observed proper- 
ties o f faint SNe II-P (e!g., SN 1997D. ICliugai fc Utrobinl 
120001 iKitaura et ail I2006D and of the low luminos- 
ity of SN 2008S-li ke transients (pos s ibly a new sub- 
class of SNe Iln, iPrieto et al.l 120081 : iThompson et al.l 
I2008D . The former ones, namely SNe 1994N, 1997D, 
1999br, 1999eu, 2001dc, and 2005cs, have been obser- 
vationally defin ed as the class of low- lum inosity Ni- 
poor SNe II-P (|Pastorello et al.l I2004L [20061 ). whose in- 
cidence is estimated to be as high as 4 — 5% of all 
SNe II. The estimated 56 Ni masses of ~ 0.002 -0.008 Mq 
for these low-luminosity SN e II-P (IZampieri et al.l 120031 : 
iPastorello et~all 12004 (20061 : iHendrv et all 120051) are in 
reasonable agreement with the present result from elec- 
tron capture supernovae. An alternative possibility is 
the origin of such supernovae from much more mas- 
sive stars (> 20 Mq) with low explosion energies, which 
suffer from fallback of freshly synthesi z ed 56 Ni (e.g., 
Turatto et all 119981: iBenetti et all 120011 : iNomoto et all 
2003t IZampieri et al. 1 120031) . This is due to the inferred 
large envelope m asses for the low-lu minosity SNe II-P 
(~ 14 - 40 M , IHendrv et all 120051 ) that favor mas- 
sive progeni tors. A recen t analy sis of the progenitors of 
SNe II-P bv lSmartt et all (|2008Q indicates, however, that 
low-luminosity supernovae with low 56 Ni production are 
likely to arise from low mass progenitors near the mini- 
mum mass limit for core-collapse supernovae. A lack of 
a-elements such as O and Mg in the case of electron cap- 



ture supernovae will be a key to sp ectroscopically distin - 
guish between these two scenarios (jKitaura et al.1 12006). 

The Crab nebula (the relic of SN 1054) is known to 
have a low kinetic energy (~ 4 x 10 49 erg. IChevalier l 
19851) and a small amount of a-elements ([Davidson et al.l 
19821 ) . An electron capture supernov a has been suggested 
to be the orig i n of the Crab nebula (|Nomoto et al.lll982l : 
iNomotd [19851: IKitaura et al.ll2006l ). Our result with the 
little production of a-elements and iron supports this 
idea. The Ni/Fe ratios (~ 1 — 2, except for the extreme 
model MX, Table 1), which are at least 20 times larger 
than the solar value (= 0.058) can be considered as an 
additional support, because they are in reasonable agree- 
ment with the reporte d high Ni/Fe ratio of the Crab neb - 
ula (~ 10 ti mes solar: lHenry||1984l: iHudgins et al.lll99(l . 

Recently, iMacAlpine fc Satterfieldl (|2008| ) investigated 
gaseous regions of the Crab nebula and inferred from 
their photoionization calculations that the abundance of 
a large component of the nebula appears to be He-rich 
and C/O > 1. According to the 8 — 10 Mq star models 
(|Nomoto et al.lll982HNomotdll984[ ). the He-rich envelope 
for M < 9.5 M Q has C/O < 1 because of the preceding 
CNO-cycle, while that for M > 9.5 Mq has C/O > 1 
owing to the 3a-reactions. The 9.6 Mq star in INomotd 
(|1984l case 2.4), for example, has carbon and oxygen 
mass fractions of 0.022 and 0.0033 in the He-burning con- 
vective layer (corres ponding solar values are 0.0025 and 
0.0066, respectively: lLoddcrs 200^). This enhanced car- 
bon abundance (ten times that of the solar value) is con- 
sistent with that reported by [MacAlpinc & Satt erfieldl 
(2008). It should be noted that stars with initial masses 
of ~ 9.5 — 10 Mq have almost identical core structures 
to our considered star of 8 .8 M ff , except for th e outer- 
most oxygen mass fraction (|Nomo tol H983 . U987h . There- 
fore, the explosion of a star with an initial mass of 
~ 9.5 — 10 M Q and with an O-Ne-Mg core can be the 
origin of the Crab remnant. 

We also note that the dredge-up of the material from 
the He-layer into the H-la yer enhances carbon in the 
envelope of the AGB star (|Nomotd 0.9871 ) . This would 
be more efficient in enhancing carbon than He-thermal 
pulses, and dust could be easily formed to induce mass 
loss. This may result in a deeply dust -enshrouded object 
such as the progenitor of SN 2008S (IPrieto et alJ l2008t 
IThompson et al.ll2008l) . For the 9.6 Mq star in iNomotol 
(|1984l case 2.4), the duration of the AGB phase is esti- 
mated to be 4 x 10 4 yr (that becomes shorter for a more 
massive case) , which is in reasonable agreement with the 
inferred dust-ensh rouded phase for SN 2 008S-like tran- 
sients (< 10 4 yr. IThompson et al.l [20081 ). This might 
also imply the mass range of the stars that end their 
lives as electron capture supernovae to be ~ 9.5 — 10 Mq, 
whose frequency, ~ 7 — 8% of all the core-collapse events, 
satisfies the constraint from our nucleosynthesis results 
(< 30%; § 4). 

6. CONCLUSIONS 

We have investigated the nucleosynthesis during the 
first 810 ms after core bounce in an explosion from a col- 
lapsing star with O-Ne-Mg core (electron capture super- 
nova) and an initial mass of 8.8 Mq. The thermodynamic 
trajectorie s are taken from the self-consistent explosion 
models of IKitaura et al.l (|2006l ). which w e re compute d 
with the initial stellar model of INomotd (|1984 1 19871) . 
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Our main conclusions can be summarized as follows. 

1. Our unmodified model (ST) results in (i) little pro- 
duction of a-elements and iron, (ii) large production of 
64 Zn, 70 Ge, and in particular, 90 Zr, and (iii) production 
of some light p-nuclei ( 74 Se, 78 Kr, 84 Sr, and 92 Mo). This 
is a consequence of the ejection of a sizable amount of 
neutron-rich matter (6 x 10~ 3 M Q with Y e = 0.46-0.49). 
If we assume this model to be representative of electron 
capture supernovae, the occurrence of this type of super- 
novae is limited to be no more than 1% of all core-collapse 
events. We do not think, however, that the production 
of 90 Zn serves as a strong constraint, because it is easily 
affected by a small variation of Y e . The z^p-process does 
not play any role for the production of p-nuclei in the 
present supernova model. 

2. The uncertainties in the nuclear EoS and the nu- 
clear reaction rates do not substantially affect the nucle- 
osynthesis results. In contrast, the effects of convection, 
which a re not include d in th e one-dimensional simula- 
tions of iKitaura et alj (|2006D . are expected to be large 
and may change the initial Y e distribution to some ex- 
tent. The overproduction of 90 Zr is moderated if the min- 
imum Y e is only 1 — 2% larger than that in the unmodified 
model ST. In this case (our model FP3) the largest over- 
production, which is observed for 64 Zn, 74 Se, and 78 Kr, 
is reduced to one tenth of that of the unmodified model. 
The robustness of the 64 Zn production against small vari- 
ations of Y e provides an upper limit to the occurrence of 
electron capture supernovae to be about 30% of all stel- 
lar core-collapse events. Electron capture supernovae can 
be significant contributors to the Galactic inventories of 
64 Zn (the most abundant isotope of Zn) and some light 
p- nuclei (e.g., 92 Mo), if the assumed slightly larger values 
of Y e were correct. 

3. The high Ni/Fe ratio (=1 — 2) and the small pro- 
duction of a-elements, as well as the low explosion en- 
ergy (1 - 2 x 10 50 erg, IKitaura et all 120061: iJanka et all 

support the hypothesis that the Crab neb- 
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ula is the remnant of an electron captu re supernova 
(|Nomoto et al.lll982t iDavidson et al.lll982f) . 

4. SN 2008S-like transients, whose progenitors are 
deeply dust-enshrouded massive stars, are likely to 
be electron capture supernovae of AGB stars. This 
implies that electron capture super novae constitute a 
newly identified sub-class of SNe Iln (|Prieto et al.ll2008t 



iThompson et al.ll2008D . 

5. The ejecta mass of 56 Ni is 0.002 - 0.004 M Q , which 
is in reasonable agreement with estimates for observed 
low-luminosity supernovae. The amount might be, how- 
ever, affected by the convective motion near the mass 
cut. Multi-dimensional studies will clarify the effect of 
convection on the production of 56 Ni as well as 90 Zr. 
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TABLE 1 

Yields in Units of Solar Masses 



Model 


Y 

1 e,mm 


56 Ni 


64 Zn 


90 Zr 




Ni/Fe 


ST 


0.464 


2.50E- 


03 


6.38E- 


-04 


1.21E- 


04 


1.65 


WH 


0.462 


4.06E- 


03 


7.31E- 


04 


1.39E- 


•04 


1.27 


RT 


0.464 


2.52E- 


03 


6.94E- 


•04 


7.83E- 


05 


1.58 


MX 


0.480 


1.67E- 


03 


1.07E- 


03 


3.32E- 


08 


3.01 


FP1 


0.468 


2.62E- 


03 


6.83E- 


04 


5.75E- 


05 


1.55 


FP2 


0.471 


2.76E- 


03 


7.08E- 


04 


1.59E- 


05 


1.46 


FP3 


0.475 


2.91E- 


03 


6.51E- 


•04 


8.04E- 


-07 


1.36 


FM1 


0.460 


2.41E- 


03 


5.83E- 


04 


1.96E- 


-01 


1.73 


FM2 


0.457 


2.32E- 


03 


5.31E- 


•04 


2.66E- 


-01 


1.82 


FM3 


0.453 


2.24E- 


03 


4.83E- 


-04 


3.11E- 


04 


1.92 



tional Wave Astronomy", and the Cluster of Excel- 
lence EXC 153 "Origin and Structure of the Universe" 
(http://www.universe-cluster.de). The computa- 
tions were done at the Rechenzentrum Garching and 
at the High Performance Computing Center Stuttgart 
(HLRS) under grant number SuperN/12758. 
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TABLE 2 

Yields of Stable Isotopes (in Units of Mq) 



Species 


ST 




FP3 




Species 


ST 




FP3 




X H . . . 


5.55E- 


05 


3.84E- 


05 


58 Ni . . 


1.79E- 


03 


2.21E- 


03 


2 H . . . 


1.98E- 


13 


1.35E- 


13 


60 Ni . . 


2.09E- 


03 


1.72E- 


03 


3 He. . . 


6.52E- 


10 


6.49E- 


10 


61 Ni . . 


5.04E- 


05 


4.62E- 


05 


4 He... 


5.12E- 


03 


5.55E- 


03 


62 Ni . . 


5.11E- 


04 


2.88E- 


01 


6 Li . . . 


6.70E- 


15 


6.57E- 


15 


64 Ni . . 


2.55E- 


07 


5.38E- 


09 


7 Li . . . 


4.52E- 


09 


5.31E- 


09 


63 Cu. . 


5.96E- 


05 


1.97E- 


05 


9 Be. . . 


2.05E- 


14 


1.99E- 


14 


65 Cu. . 


1.56E- 


05 


4.59E- 


06 


10 B . . 


2.72E- 


14 


2.45E- 


14 


64 Zn. . 


6.38E- 


04 


6.51E- 


01 


X1 B . . 


3.02E- 


09 


3.34E- 


09 


66 Zn. . 


4.54E- 


04 


2.16E- 


05 


12 C . . 


1.18E- 


06 


1.48E- 


06 


67 Zn. . 


4.44E- 


06 


2.03E- 


06 


13 c . . 


6.82E- 


09 


7.09E- 


09 


68 Zn. . 


2.89E- 


05 


3.90E- 


05 


14 N . . 


5.54E- 


08 


4.08E- 


08 


70 Zn. . 


9.72E- 


12 


2.14E- 


14 


15 N . . 


1.02E- 


07 


8.25E- 


08 


69 Ga. . 


1.93E- 


-06 


5.68E- 


07 


16 o . . 


1.03E- 


07 


1.13E- 


07 


71 Ga. . 


3.96E- 


07 


7.64E- 


08 


17 o . . 


5.17E- 


11 


5.33E- 


11 


70 Ge. . 


8.29E- 


05 


6.61E- 


06 


18 o . . 


7.16E- 


09 


5.41E- 


09 


72 Ge. . 


5.05E- 


-06 


7.10E- 


07 




1.02E- 


08 


8.95E- 


09 


73 Ge. . 


4.11E- 


07 


4.32E- 


08 


20 Ne. . 


3.12E- 


08 


4.08E- 


08 


74 Ge. . 


5.16E- 


09 


6.86E- 


12 


91 IV T 


1.63E- 


10 


2.00E- 


10 


7b Ge. . 


2.91E- 


-14 


4.90E- 


20 


22 Ne. . 


5.37E- 


09 


5.82E- 


09 


75 As . . 


3.35E- 


07 


1.84E- 


08 


23 Na. . 


3.33E- 


10 


4.46E- 


10 


74 Se . . 


5.13E- 


06 


9.28E- 


07 


24 Mg . 


1.80E- 


08 


2.53E- 


08 


76 Se . . 


9.05E- 


-06 


3.39E- 


08 


25 Mg . 


5.30E- 


09 


1.05E- 


08 


77 Se . . 


2.77E- 


07 


7.99E- 


09 


26 Mg . 


4.40E- 


08 


4.30E- 


08 


78 Se . . 


1.32E- 


07 


5.41E- 


11 


27 A1 . . 


3.03E- 


09 


3.32E- 


09 


80 Se . . 


1.49E- 


10 


2.87E- 


11 


28 Si. . . 


5.44E- 


08 


8.45E- 


08 


82 Se . . 


8.40E- 


16 


0.00E+00 


29 Si. . . 


1.11E- 


08 


1.66E- 


08 


79 Br . . 


2.08E- 


07 


1.15E- 


08 


30 Si. . . 


6.82E- 


08 


6.13E- 


08 


81 Br . . 


1.23E- 


07 


3.09E- 


09 


31 p 


2.10E- 


08 


2.82E- 


08 


78 Kr . . 


7.06E- 


07 


2.58E- 


07 


32 s ... 


1.98E- 


07 


2.76E- 


07 


80 Kr . . 


2.88E- 


06 


4.71E- 


08 


33 s ... 


3.37E- 


08 


3.89E- 


08 


82 Kr . . 


1.08E- 


06 


2.53E- 


08 


34 s ... 


2.27E- 


07 


2.20E- 


07 


83 Kr . . 


2.03E- 


07 


5.58E- 


09 


36 s ... 


4.09E- 


12 


5.03E- 


12 


84 Kr . . 


3.57E- 


08 


6.87E- 


12 


35 C1 . . 


1.78E- 


07 


2.03E- 


07 


86 Kr . . 


1.24E- 


-11 


3.48E- 


18 


37 C1 . . 


1.05E- 


07 


1.14E- 


07 


85 Rb. . 


1.40E- 


07 


4.29E- 


09 


36 Ar . . 


1.11E- 


06 


1.41E- 


06 


87 Rb. . 


4.21E- 


09 


1.39E- 


15 


38 Ar . . 


2.62E- 


07 


3.01E- 


07 


84 Sr . . 


3.10E- 


07 


5.87E- 


08 


40 Ar . . 


9.90E- 


11 


1.24E- 


10 


86 Sr . . 


1.29E- 


06 


6.94E- 


09 


39 K . . 


4.07E- 


07 


4.83E- 


07 


87 Sr . . 


2.92E- 


07 


2.72E- 


09 


40 K . . 


4.98E- 


10 


5.98E- 


10 


88 Sr . . 


2.74E- 


06 


3.24E- 


08 


41 K . . 


8.14E- 


08 


8.53E- 


08 


89 y 


3.92E- 


06 


1.25E- 


08 


40 Ca. . 


3.55E- 


06 


3.93E- 


06 


90 Zr . . 


1.21E- 


04 


8.04E- 


07 


42 Ca. . 


5.23E- 


07 


6.98E- 


07 


91 Zr . . 


9.96E- 


07 


6.61E- 


08 


43 Ca. . 


1.77E- 


07 


2.27E- 


07 


92 Zr . . 


7.81E- 


09 


1.35E- 


10 


44 Ca. . 


5.96E- 


06 


7.17E- 


06 


94 Zr . . 


3.89E- 


16 


1.97E- 


20 


46 Ca. . 


2.44E- 


1.5 


1.25E- 


15 


96 Zr . . 


0.00E+00 


0.00E+00 


48 Ca. . 


0.00E+00 


0.00E4 


-00 


93 Nb. . 


2.05E- 


08 


2.07E- 


09 


45 Sc . . 


1.41E- 


07 


1.46E- 


07 


92 Mo . 


7.02E- 


07 


2.73E- 


07 


46 Ti . . 


7.18E- 


07 


8.29E- 


07 


94 Mo . 


1.42E- 


08 


2.32E- 


10 


47 Ti . . 


9.24E- 


07 


1.17E- 


06 


95 Mo . 


7.96E- 


11 


1.16E- 


11 


48 Ti . . 


8.36E- 


06 


1.01E- 


05 


96 Mo . 


3.94E- 


12 


6.44E- 


11 


49 Ti . . 


7.47E- 


07 


7.64E- 


07 


97 Mo . 


1.46E- 


12 


2.09E- 


13 


50 Ti . . 


3.35E- 


12 


7.44E- 


14 


98 Mo . 


6.96E- 


18 


0.00E+00 


50 V . . 


1.94E- 


09 


2.97E- 


10 


100 Mo 


0.00E+00 


0.00E+00 


bi V . . 


1.95E- 


06 


2.23E- 


06 


Qfi-r-. 

Jt> Ru. . 


2.59E- 


-11 


1.58E- 


11 


50 Cr . . 


1.41E- 


06 


1.64E- 


06 


98 Ru.. 


2.86E- 


-12 


6.57E- 


11 


52 Cr.. 


1.50E- 


05 


1.82E- 


05 


"Ru.. 


1.70E- 


-14 


8.73E- 


16 


53 Cr.. 


1.04E- 


06 


1.22E- 


06 


100 Ru. 


7.23E- 


15 


4.99E- 


16 


54 Cr . . 


1.26E- 


08 


1.86E- 


09 


101 Ru. 


2.58E- 


16 


2.17E- 


17 


55 Mn . 


3.41E- 


06 


2.88E- 


06 


102 Ru. 


9.02E- 


19 


0.00E+00 


54 Fe . . 


3.22E- 


06 


3.64E- 


06 


104 Ru. 


0.00E+00 


0.00E+00 


56 Fe . . 


2.52E- 


03 


2.92E- 


03 


103 Rh. 


4.39E- 


18 


1.01E- 


19 


57 Fe . . 


1.80E- 


04 


2.13E- 


04 


102 Pd. 


3.81E- 


16 


1.65E- 


17 


58 Fe 


7.25E- 


08 


2.18E- 


08 


104 Pd. 


3.19E- 


18 


4.54E- 


21 


59 Co'.'. 


8.61E- 


05 


1.08E- 


04 


105 Pd. 


1.52E- 


-20 


0.00E4 


-00 



TABLE 3 

Yields of Unstable Isotopes (in Units of Mq) 



Species 


ST 




FP3 


Species 


ST 




FP3 


22 Na. . 


5.37E- 


09 


5.81E-09 


60 Fe . . 


1.40E- 


14 


4.79E-16 


26 Al . . 


3.27E- 


08 


2.97E-08 


56 Ni . . 


2.50E- 


03 


2.91E-03 


41 Ca.. 


8.05E- 


08 


8.43E-08 


57 Ni . . 


1.77E- 


04 


2.11E-04 


44 Ti . . 


5.96E- 


06 


7.17E-06 


92 Nb.. 


6.35E- 


09 


1.34E-10 
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